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Abstract 

CFD calculations using high-performance parallel comput- 
ing were conducted to simulate the pre-stall flow of a tran- 
sonic compressor stage, NASA compressor Stage 35. The 
simulations were run with a full-annulus grid that models the 
3D, viscous, unsteady blade row interaction without the need 
for an artificial inlet distortion to induce stall. The simulation 
demonstrates the development of the rotating stall from the 
growth of instabilities. Pressure-rise performance and pressure 
traces are compared with published experimental data before 
the study of flow evolution prior to the rotating stall. Spatial 
FFT analysis of the flow indicates a rotating long-length dis- 
turbance of one rotor circumference, which is followed by a 
spike-type breakdown. The analysis also links the long-length 
wave disturbance with the initiation of the spike inception. 
The spike instabilities occur when the trajectory of the tip 
clearance flow becomes perpendicular to the axial direction. 
When approaching stall, the passage shock changes from a 
single oblique shock to a dual- shock, which distorts the per- 
pendicular trajectory of the tip clearance vortex but shows no 
evidence of flow separation that may contribute to stall. 

Introduction 

Advances in gas turbine engine technology have made pos- 
sible the production of highly efficient engines. However, the 
useful operating range of a gas turbine engine is limited by the 
stable operation of the compressor. The limiting instability can 
be characterized by either stall or surge. Surge is the zeroth 
order planar oscillation mode resulting from overall imbalance 
of the entire propulsion system, thus, analysis would have to 
include full volume and ducting. This violent system instabil- 
ity is always preceded by a rotating stall. Surge and stall can 
cause destructive damage to engines and should be avoided. 
Understanding the physical process of the onset of rotating 
stall can aid the design of better control method to stabilize the 
system and enhance engine performance. 

Understanding stall inception has been an active research 
effort for the past five decades. The events leading to rotating 
stall have been traditionally classified according to two differ- 
ent types of wave disturbances moving around the annulus: 


long-length (modal) and short-length (spike) waves. Long- 
length wave disturbances are those with length on the order of 
the rotor circumference. This type of disturbance is caused by 
the characteristic frequencies of the compression system. 
Theories to model the long-length instability were developed 
by Moore (1984), and Moore and Greitzer (1986). The as- 
sumption is that the disturbance length scale is much larger 
than one blade pitch, and that the entire blade row can be 
modeled, both spatially and temporally, without the knowl- 
edge of the detailed flow inside the blade passages. This ap- 
proach predicts the existence of long-length waves leading to 
stall. McDougall (1988, 1989) demonstrated the existence of 
the modal waves through experiments, and Gamier et al. 
(1991) supported the validity of the Moore and Greitzer model 
with detailed measurements. Unsteady pressure or axial veloc- 
ity measurements from probes placed around the annulus are 
usually used to track the development of modal waves. Modal 
oscillations often have an axial extent from inlet to exit. It is 
typically a 2D (0,x) phenomenon, which is not an early form 
of stall cell but instead represents harmonic oscillations of the 
flow field. It usually appears when the shape of the total-to- 
static pressure rise characteristic becomes horizontal. 

The second type of stall inception mechanism is the short- 
length, or spike disturbances. These are disturbances with 
length scale on the order of blade passage breadth. Spike can 
be viewed as embryonic stall cells with flow breakdown in 
local regions (Camp & Day, 1998), it is a 3D phenomenon 
whose development depends on the flow stmcture within the 
blade passage. A common scenario for the spike stall incep- 
tion process is that localized, 3D pockets of instability in mul- 
tiple blade passages empt first to form small sized multi-cell 
rotating stall. These stall cells rotates at 30 to 85% rotor speed 
initially and then merge into a single rotating stall cell within a 
few rotor revolutions (Day, 1993). The final speed of the con- 
solidated stall cell is nearly half or less of the rotor speed 
(Cumpsty, 2004). A spike disturbance mostly occurs when the 
slope of the pressure rise characteristic is still negative. That 
is, the pressure still rises as the mass flow reduces. 

The above discussions indicate that the slope of the pres- 
sure rise characteristic may correlate with the stall inception 
type. Spike and mode can also be the different aspects of the 
stall inception process and can exist in the same machine 
(Day, 1993). Day indicated that modal waves may be measur- 
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able in the machine for as many as 200 rotor revolutions be- 
fore stall. A spike disturbance, on the other hand, once 
formed, can lead to rotating stall quickly, usually within a few 
rotor revolutions. 

Numerical simulation using CFD techniques has been 
shown as an alternative way to analyze aerodynamic problems 
and is a complement to experiments in the understanding of 
complex flow phenomena. Flow structure can be examined in 
detail where experimental measurements are difficult or im- 
possible. During the last ten years steady-flow simulations 
have provided an increasingly accurate prediction of the flow 
up to the point of compressor stall. Attempts to study stall 
through unsteady simulations of a periodic subset of the blades 
in a compressor blade row (Hoying et al., 1988, Vo et al., 
2005, Davis and Yao, 2006) or through reduced-order un- 
steady flow models (Gong et al. 1999) have reached valuable 
findings of the flow structure leading to stall. However, since 
the temporal flow field variations that occur during stall incep- 
tion are not harmonics of blade passing frequency, a full- 
annulus model should be used to remove the constraint of the 
periodic model in order to study the transition from a steady 
flow into the unsteady stalling flow, and to resolve flow fea- 
tures of length scale on the order of the rotor circumference. 
Such simulations have been done two-dimensionally (He, 
1997, Saxer-Felici et al., 2000), but stall and the instabilities 
leading to stall are inherently three-dimensional phenomena. 
Predicting stall phenomena may thus require three- 
dimensional unsteady full-annulus simulations, which is a 
daunting computational endeavor requiring considerable com- 
putational resources. Due to the rapid advance of high per- 
formance computing in the past decade, advanced CFD codes 
have benefited from the almost unlimited computing potential 
of parallel computing. Recently, Chen et al. (2006) demon- 
strated the ability to conduct a full-annulus 3D unsteady com- 
putation of an entire compressor stage while Hah et al. (2006) 
also conducted a full-annulus 3D computation of an isolated 
rotor. This trend indicates that large-scale simulation using 
parallel computing is within reach. 

Previous comparisons of the time-averages of unsteady 
phase-lag and partial annulus predictions to measurements 
have shown the validity of the TURBO code to capture the 
ID, and 2D characteristics of this compressor (Hathaway et 
al., 2004). The intent of this work is to use the available ex- 
perimental data to provide assessment of the code and to show 
evidence that the computational model reasonably captures the 
inherent behavior of compressors as they are throttled into 
stall. It is recognized that some features of the actual tested 
compressor are not captured by the computational model, for 
example the leakage gap at the rotor/stator hub interface, the 
actual blade tip profile, and the geometry of the upstream and 
downstream ducting between the plenum and throttle valve 
locations. Although these features might influence the actual 
compressor stall, they are not considered requisite to the un- 
derstanding of the fundamental stall inception characteristics. 


Axial Compressor, Stage 35 

NASA research single-stage compressor Stage 35, repre- 
sentative of the advanced transonic core compressor, is used in 
this investigation. A schematic of the compressor test rig and 
compressor geometry is provided in figure 1 and the design 
parameters are provided in Table 1. Stage 35 produces 1.8 
total pressure ratio at a mass flow rate of 20.2 kg/sec at the 
design speed of 17189 rpm. Details of the Stage 35 geometry, 
operating conditions, and the experimental results are pro- 
vided by Reid and Moore (1978). 


TABLE 1.— DESIGN PARAMETERS OF STAGE 35 


Rotor rpm at 100% Speed 17188.7 

Tip Speed (m/s) 454.456 

Hub/Tip Radius 0.7 

Rotor Aspect Ratio 1.19 

Stator Aspect Ratio 1 .26 

Number of Rotor Blades 36 

Number of Stator Blades 46 



Figure 1. — NASA high-speed compressor Stage 35. 
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Relevant Experimental Work 

Time-resolved pressure measurement of Stage 35 near stall 
was reported by Bright et al. (1999). The experiment was 
performed at 85% design speed. Several configurations, in- 
cluding a clean-inlet case, were reported in that paper. For the 
interest of this study we will focus on the clean-inlet case. The 
compressor was transitioned into stall using a continuous 
throttle maneuver. Time-resolved pressure was recorded by 
high-frequency-response transducers at eight locations equally 
spaced around the annulus. Figure 2 reproduces the pressure 
traces of the clean-inlet case for 1 80 revolutions prior to stall 
reported in that paper (the original data is not available to us). 
Based on the pressure disturbance, Bright et al. commented 
that the compressor with the clean inlet stalls due to modal 
instabilities, and at 40 revolutions before stall there appear to 
be spikes on top of the modal wave. They concluded that 
“This indicates that our high-speed stage under normal stall- 
ing condition is considered a ‘modal' machine, since modal 
instabilities dominate the pre-stall behavior with some added 
pip instabilities present just before stall ” 

In the following we take a closer look at this experimental 
data. The shapes of the modal disturbances indicated in 
figure 2 resemble sinusoidal waves with a temporal period of 
10 rotor revolutions. The speed of the modal disturbances, as 
indicated by the slope of line A, connecting the peaks of the 
pressure wave, is zero. This indicates the modal disturbance is 
caused by axisymmetric longitudinal system disturbances. The 
system includes all ducting and volumes of the experimental 
facility. This modal disturbance is the surge mode that oscil- 
lates at a frequency of 10% rotor speed — equivalent to the 10 
rotor revolutions period it covers. Unlike a typical rotating 
modal wave, the surge mode does not rotate. The surge mode 
is a compression system wide resonance that cannot be mod- 
eled with the current simulation, because the computational 
grid does not include the ducting of the test facility, which is 
quite substantial, as shown in figure 1. The surge mode oscil- 
lates over 10 rotor revolutions and is characterized by a length 
scale of 10 rotor circumferences, or, more than one rotor cir- 
cumference. 

As the compressor throttles closer to stall at about 1 0 revo- 
lutions before stall, there is a secondary disturbance, much like 
a spike in the time scale of figure 2. Looking more closely, 
this secondary disturbance covers a circumferential distance of 
one rotor circumference, as indicated in the enlarged image. It 
sits on top of the surge mode and rotates around the annulus at 
about 80% rotor speed, indicated by line B. Because the length 
scale of this disturbance is the entire annulus, which is much 
larger than that of a typical spike, it is classified as a mode. 
For discussions in the rest of the paper, this rotating long- 
length disturbance is also referred to as a rotating modal dis- 
turbance. This long-length disturbance is followed by the 
eventual flow breakdown. Due to the picture quality, we are 
unable to discern the transition between the start of the rotat- 
ing long-length disturbance and the final breakdown. The 
long-length disturbance has a length scale of one rotor 



Rotor Revolutions 

Figure 2. — Experiment pressure traces for clean inlet, 
85% design speed, Stage 35, Michelle Bright et al. (1999). 



Figure 3. — Comparison of characteristics at 
85 and 100% speed, Stage 35. 


circumference, or, less than one rotor circumference. This is in 
contrast to the more than one rotor circumference of the sur- 
gemode. The rotating long-length disturbance, because of its 
length scale, can only be resolved by a full annulus simulation. 

The simulation reported in this work is part of an overall 
| effort to examine steady tip injection stall flow control for 
Stage 35 at design speed with the clean inlet. As such, the 
simulation was conducted at design speed. To our knowledge, 
there is no pressure trace data available at design speed as that 
shown in figure 2. Therefore, we will use the 85% speed case 
to compare with the simulated results. While the speeds are 
different, the two cases are both transonic compressors and 
share common features as we explore the flow behavior prior 
to stall. The experimental characteristics at 85 and 100% 
speed are compared in figure 3. Both case stalls at the peak of 
static pressure rise. More importantly, their slopes near stall 
are similar. These similarities justify the use of 85% speed for 
comparison. 
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Numerical Methodology 

TURBO is a physics-based simulation tool for multistage 
turbomachinery. The solver computes the fluid conservation 
laws without ad hoc modeling of any flow phenomena other 
than models required for turbulence. This code solves the 
unsteady Reynolds-averaged Navier-Stokes equations and a 
decoupled k-s turbulence model. The code is implemented in a 
portable, scalable form for distributed-memory parallel com- 
puters using MPI message passing. The parallel implementa- 
tion employs domain decomposition and supports general 
multi-block grids with arbitrary grid-block connectivity. The 
solution algorithm is a Newton iterative implicit time-accurate 
scheme with characteristics-based finite-volume spatial discre- 
tization. The Newton subiterations are solved using a concur- 
rent block-Jacobi symmetric Gauss-Seidel (BJ-SGS) 
relaxation scheme. Because all of the fundamental fluid me- 
chanics are computed, the code is capable of capturing the 
nonlinear characteristics of the flow fields of interest. With the 
actual modeling of blade rows in relative motion, this code is 
capable of computing the unsteady interactions between blade 
rows. Details of the flow solver can be found in Chen and 
Whitfield (1993). The approach to parallelization for large- 
scale, complex problems is discussed in Chen and Briley 
(2001). In the following we discuss several issues applied to 
the code in the simulation of rotating stall. 


lows variation of exit static pressure to match the compressor 
exit mass corrected to the exit total condition. 

Tip Clearance Model 

The tip gap was modeled by the tip clearance model of Kir- 
tley (1990). This model ensures the conservation of mass and 
momentum through the tip gap without the need of a tip clear- 
ance grid. A gridding strategy suggested by VanZante et al. 
(2000) is adopted in this study where 10 points are placed in 
the tip gap region. Although the detailed physics of the devel- 
opment of the leakage vortex are not captured the tip clearance 
model was considered adequate for predicting the first order 
effects of the leakage vortex strength, extent, and direction 
(Chima, 1998). 

Gridding 

This compressor stage consists of 36 rotor blades and 46 
stator blades. For the purpose of conducting simulations with 
tip injectors in another effort, we built a three-blade-row grid 
as shown in figure 4. The injector row has twelve injectors 
(not activated in this study) that are equally spaced in the 
circumferential direction. The grid size of the rotor is 151(x), 
71(r), and 56(0) with 81(x) and 61(r) points on the blade. The 
grid size of the stator is 141(x), 71(r), and 79(0) with 81(x) 
and 53(r) points on the blade. The rotor tip clearance is 1% of 


Full Compressor Stage Simulation 

The intent to simulate the inception and maturation of long- 
wavelength circumferential modes, with wave numbers as low 
as unity, requires full annulus simulations. A time-accurate 
sliding interface (Chen and Barter, 1998) is implemented, in 
which pressure waves in all directions (axial, radial, and 
circumferential) are passed unaltered between the rotor and 
the stator. This attribute is necessary to examine the temporal 
and spatial responses associated with blade row interaction 
during stall inception. 

Inlet Boundary Condition 

An isentropic inlet condition (Chen and Whitfield, 1993) is 
applied, where the span-wise total condition is prescribed. 
This boundary condition preserves the upstream total condi- 
tions at test facility level, but can potentially reflect out-going 
pressure waves back to the computational domain. 

Exit Boundary Condition 

The traditional choice of exit boundary for stable flow is the 
radial equilibrium condition. This boundary condition uses a 
preset exit static pressure. In near-stall or during-stall condi- 
tions, compressor pressure rise drops and a preset exit pressure 
cannot match the pressure drop. In this study we relieve the 
situation by the use of a “choked” throttle model that specifies 
corrected mass flow at the exit. This boundary condition al- 



Figure 4. — Three-blade-row Grid Model 
for the Stage 35 simulation. 
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rotor tip chord. The grids are replicated around the annulus to 
generate a full-annulus grid with a total of 67 million grid 
points. The grid is then partitioned into 328 blocks. The simu- 
lations were run on an IBM P4 cluster at NAVO MSRC using 
328 CPUs with 24 hr (wall clock) computation time for one 
rotor revolution. In the simulation of the development of rotat- 
ing stall, 10 revolutions were needed. Ideally, each operating 
point should be run at least 10 revolutions to ensure conver- 
gence, especially for points close to stall. The time for each 
speedline would then be 0.5 millions CPU-hours (6 operating 
points). 

Constrained by the computational resources for the full an- 
nulus simulations, very short upstream and downstream com- 
putational domains were used here. This choice prompts two 
caveats: 1) The exclusion of the test facility volume effect 
prevents the modeling of the surge mode. Any modal activity 
shows up in the simulation should be a result of the compound 
effect of the blades (both rotor and stator if the interaction of 
the two matters) and the limited volume effect present in the 
ducting of the current grid. 2) Limiting the axial extent of the 
computational domain can potentially disrupt long length 
modes due to the reflective inlet boundary. Because pressure 
waves (with length X) decay e _ZA through a distance L, the 
common practice is to locate inlet boundary at least one rotor 
circumference away from the rotor so that pressure wave of 
one rotor circumference decays sufficiently before reaching 
the inlet. Typically, the modal instability shows up hundreds 
of revolutions prior to stall. Complete simulation of the evolu- 
tion of such modal instability is beyond the computational 
resources available in this study. Our intent, therefore, is to 
examine the development of the modal disturbances as they 
first emerge before reaching the inlet, rather than the complete 
resolution of the modal stall process. In this aspect, a shorter 
upstream domain is justified. 

Results and Discussions 

The difference of the throttling processes in the simulation 
and the experiment is that in the simulation the exit boundary 
condition is changed in discrete steps, while in the experiment 
the compressor transitioned to stall through a continuous throt- 
tling maneuver. In this simulation, a new throttle setting re- 
starts from the solution of a previous throttle setting. Due to 
computational resources, each throttle setting was only run for 
a few revolutions before moving to the next. Generally, it was 
sufficient to reach convergence within 4 revolutions while 
away from stall. As stall was approached, more (6 or more) 
revolutions were required to assure convergence. 

Stage Characteristics 

The computed and experimental characteristics (speedlines) 
are shown in figure 5. The data used in this figure is the ratio 
of exit static pressure to inlet total pressure versus mass rate 
corrected to inlet total condition. The computed speedline is 


generated by incrementally closing the throttle setting (de- 
creasing exit corrected mass flow) from A to F. The overall 
shape of the computed characteristic is similar to the experi- 
mental data, so is the slope of the characteristic near the peak 
pressure (throttle point D). Since the slope of the compressor 
characteristic correlates with the type of stall inception, we 
consider the simulation reasonable to capture the stall incep- 
tion mechanisms. 

Point D in figure 5 is the last stable throttle setting. Points E 
and F both restart from Point D and develop into rotating stall. 
In both cases the exit pressure changes from the peak at D into 
a transitional peak and then drops to a much lower level dur- 
ing rotating stall, as shown in figure 6(a). This demonstrates 
the throttling effect of the corrected mass exit boundary condi- 
tion by allowing the exit pressure to vary as the flow condition 
changes. The pressure ratios shown in figure 5 for E and F are 
the transitional peak pressures before entering rotating stall. 
Point D* is the approximate stall point from the experiment. 
The actual simulated stall point should occur between D and 
E. The actual stall point depends on many issues, such as tip 
clearance height, deviation of real blade geometry from de- 
signed geometry, as well as the type of turbulence modeling in 
the numerical simulation. These details may not be correctly 
represented in this simulation. Even though the exact stall 
point is important, for the intent of this paper the simulation is 
considered adequate to examine the general physical trends 
leading to stall. 

In the following sections, we first present the entire evolu- 
tion of the rotating stall as the throttle setting changes from D 
to F. This is then followed by a closer examination of the pre- 
stall flow development from D to E, with E being closer to the 
stable point D. 



Figure 5. — Stage 35 computed and 
measured speedline, design speed. 
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revolution 
(a) Pressure 



Figure 6. — Pressure and mass variation during 
rotating stall at throttle setting F. 



Figure 7. — Time history of static pressure variation at 8 loca- 
tions around the annulus located 44% chord ahead of 
the rotor of Stage 35, throttling setting F. 


Evolution of the Rotating Stall 

This result was reported in Chen et al. (2006) and is out- 
lined here for completeness. The rotating stall was captured by 
the time variation of static pressure at eight pressure transduc- 
ers circumferentially spaced upstream of the rotor. The probes 
are located 44% chord upstream of the rotor leading edge at 
98% span. The circumferential locations are 10°, 70°, 100°, 
160°, 190°, 250°, 280°, and 340°. The history of pressure 
variations are shown in figure 7, in which the pressure levels 
are offset by the probe circumferential locations. Because the 
volume effect is not modeled in the simulation, the pressure 
disturbance of figure 7 should be viewed as a filtered version 
of figure 2 by removing the surge mode oscillations. 

Small, periodic pressure variations caused by the rotor 
blade passing, as indicated by the thickness of the curves, are 
present in all curves of figure 7. A smooth disturbance is visi- 
ble at t=0.7T (T = time for one rotor revolution) traveling at 
full rotor speed. This disturbance, analyzed in further detail 
later, has a length scale of one rotor circumference, thus is 
classified as a modal disturbance. The rotating modal distur- 
bance continues to grow in magnitude but the speed slows to 
84% rotor speed in approximately one revolution. After 2T, it 
transforms into several spike disturbances around the annulus, 
each with a smaller length scale of a few blade passages. The 
spike disturbances then grew in magnitude and slowed down 
to 43% in 4 revolutions. Inlet and exit mass flow histories 
during the same period are shown in figure 6(b), in which the 
mass flow remained relatively flat in the first two revolutions 
at 18 kg/s before it dropped sharply and underwent large 
variations of 14% of the peak mass flow. 

In figures 8 and 9, the entropy on a surface of revolution 
near the casing is shown at two instants in time. The full- 
annulus surface is illustrated by two half- annulus surfaces. At 
t=3.8T, figure 8, the early phase of the stall inception, three 
separate high entropy (red color) regions can be seen, with the 
axial extent covering from one chord upstream to aft of the 
stator leading edges. The radial extent of the stall cell is shown 
in the axial cut-plane at 16 % chord before the rotor. Again, 
the high entropy region shows the stall cell extending about 
30% span from the casing. At 1=1 0T, figure 9, the stall cell 
evolves into the final stage as its rotating speed became stabi- 
lized. The three-cell stall has merged into a single-cell stall 
with the axial extent covering all the way from the inlet of the 
computation domain to aft of the stator trailing edge. The 
radial extent of the stall cell has now extended about 70% span 
from casing. 

The evolution of the simulated rotating stall cell is similar 
to the spike stall inception mechanism commonly observed in 
that localized, three-dimensional pockets of instability in mul- 
tiple rotor passages break up to form the initial phase of a 
rotating stall cell, usually with several small-sized stall cells. 
These small stall cells move fast at the beginning and then 
merge into a single rotating stall cell within a few revolutions 
before slowing down to about half of the rotor speed as the 
rotating stall becomes fully developed. 


NAS A/TM— 2008-2 15163 6 






16% chord upstream of rotor leading edge 


Figure 8. — Entropy of Stage 35 during 
stall inception process, t=3.8T. 



Left half annulus Right half annulus 
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16% chord upstream of rotor leading edge 

Figure 9. — Entropy of Stage 35 during 
stall inception process, t=10T. 


In the previous section we gave a brief description of the 
overall evolution of the rotating stall using throttle setting F. 
In the following sections we will examine the process that 
leads to the initial flow breakdown, i.e., the stall inception 
process. Restarting from the stable solution of D, the simula- 
tion was repeated with throttle setting E, the mass flow history 
during the transition from D to E is shown in figure 10. The 
simulation shows that after four revolutions after the restart 
the mass flow starts to decrease, leading to the rotating stall. 
Three time stamps marked El (1.7T), E2 (3.2T), and E3 
(4.3T) will be referenced in the following discussions. 

Development of Long-Length Disturbance 

Both experimental (Spakovszky et al., 1999 and Suder et al. 
2001) and computational studies (Chen et al., 2006, Davis and 
Yao et al., 2006) indicate the rotating stall of this compressor 
initiates in the rotor tip region, so the analysis will focus on 
the rotor tip. Instantaneous pressure profiles around the annu- 
lus near the rotor tip leading edge are shown in figure 11. 


The three curves correspond to different phases at throttle 
setting D and E during the transition from stable operation to 
stall (fig. 10). Phase 1 is the last stable point before stall, 
which will be used as the baseline for comparison. Phases 2, 3 
are the same throttle setting E, but at two different instants: 
one before stall (time El), the other at the beginning of stall 
(time E2). Spatial FFTs of the pressure profile are shown in 
figure 12. The most dominant harmonic is the 36 th mode (its 
magnitude is at least one order of magnitude larger than oth- 
ers), which corresponds to the number of rotor blades. This 
disturbance is caused by the periodic flow features between 
blades and thus has the length scale associated with the rotor 
blade pitch. Not shown here is the strength of the 46 th mode, 
which is negligible. This indicates the influence of the stator 
on the potential field at rotor leading edge can be ignored. 
Therefore we will focus the study only on the rotor. 

Modal disturbances, if any, should exist in the low harmon- 
ics with a longer length scale that is on the order of the rotor 
circumference. At phase 1 (D), only the 36 th mode prevails. 
This indicates that when away from stall, the only disturbance 
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magnitude (log scale) normalized static pressure 



Figure 10. — Mass flow history at throttle setting E. 



Figure 11. — Instantaneous pressure profile at rotor tip. 



Figure 12. — FFT of figure 11. 


to the otherwise uniform annulus flow is that from the periodic 
rotor passing and there is no modal activity. At phase 2 (time 
El), the early phase of the rotating stall, there exists an identi- 
fiable 1 st mode along with the 36 th mode. There are also 
weaker but identifiable modes between the 1 st and the 36 th 
modes. As the simulation proceeds, bringing the compressor 
even closer to stall at phase 3 (time E2), the magnitude of the 
1 st mode increases even more. The magnitudes of the modes in 
between remain relatively the same as those in the previous 
two phases, but fluctuate over the spectrum. The only mode 
that grows persistently during the pre-stall phase is the 1 st 
mode. This disturbance rotates at 100% rotor speed. This 
observation is the same as that of throttle setting F, where the 
speed of rotating mode is indicated in the pressure history of 
figure 7. This finding is similar to the experiment (see fig. 2) 
in that the less than one rotor circumference disturbance, 
which sits on top of the surge mode (not modeled in this simu- 
lation), starts to rotate at 80% rotor speed right before stall. It 
suggests perhaps that the simulation is adequately modeling 
the less than one rotor circumference modal disturbance be- 
fore the compressor enters rotating stall. 

A possible explanation of the origin of the 1 st mode can be 
derived from the following observations: 1) while the grid in 
our simulation does not include the complete ducting system, 
it does include the limited ducting that contains the rotor and 
the stator, 2) the speed of the rotating mode is 100% rotor 
speed; in another words, the 1 st mode is fixed to the rotor. This 
makes it likely that the disturbance is associated with the rotor. 
The two observations suggest that the 1 st mode may originate 
from a compound effect of the rotor blades and the limited 
volume effect present in the ducting of the current grid. 

Development of Mode and Spike 

To better understand the development of the 1 st mode and 
the spike disturbance, we focus on the evolution of the flow 
during the transient right before stall. Since stall is likely to be 
associated with the reverse flow (negative axial flow) we will 
monitor the development of the negative axial velocity of the 
flow during the transition. Figure 13(a) to (c) illustrates the 
axial velocity of the stage at three radial cut-planes at the three 
time instants El, E2, E3 during the transient. The three cut- 
planes are located at: 4 clearance heights (4s) below casing 
(cut-plane 1), 1.5 clearance heights ((1.5s) 1. below casing 
(cut-planes 2), and 0.5 clearance heights (0.5s) below casing 
(cut-plane 3). Figure 13(d) illustrates the front views of cut- 
plane 3 at the three time instants. The contours of the axial 
velocity have been truncated so that only the negative velocity 
has a color. All positive velocities are colored white. 

At time El, figure 13(a), about 2.5 revolutions before 
breakdown, the reverse flow is seen in all rotor passages with 
a radial coverage of at least 4 clearance heights below the 
casing, as indicated by the reverse flow being visible on all 
cut-planes. The rotor reverse flow is formed by the tip clear- 
ance vortex flow, which will be explained in detail later. The 
boundary and magnitude of the reverse flow region depicts the 
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(a) Time El (2.5 T before rotating stall) 



(b) Time E2 (beginning of rotating stall) 



(c) Time E3 (during rotating stall) 

Cut-plane 1 (4s) Cut-plane 2 (1.5s) Cut-plane 3 (0.5s) 

Figure 13. — (a), (b), and (c). Negative axial velocity at 3 cut-planes, 3 time instants. 
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Time El Time E2 Time E3 

Figure 13. — Concluded, (d) Front views of cut-plane 3 (0.5s) at three time instants, negative axial velocity. 


boundary and the strength of the tip vortex. The patterns of the 
reverse flow are not the same in every passage, with the re- 
verse flow zone larger on the back side of the view (210°). 
This can be readily seen on cut-planes 1 and 2. The overall 
shape of rotor reverse flow corresponds well with the harmon- 
ics of figure 12 in that the 1 st mode co-exists with the 36 th 
mode. The speed of the 1 st mode is the same as the 36 th mode, 
both rotate at the rotor speed. The flow at this phase has al- 
ready shown the instability. 

At time E2, by comparing the front view of cut-plane 3 at 
time El and E2 (figure 13(d)), we see the front boundary of 
the reverse flow moves forward and becomes perpendicular to 
the axial direction. The magnitude of the negative velocity at 
time E2 is also higher than that at time El as is shown by 
more green spots at time E2. This is consistent with the short- 
length stall inception mechanism described by Hoying et al. 
that spike stall initiates when the tip clearance vortex trajecto- 
ries align perpendicularly to the axial direction. The larger 
reverse flow zone now rotates to the front side of the view 
(0°), figure 13(b), better seen on cut-planes 1 and 2. Also seen 
on the two cut-planes is the appearance of reverse flow in the 
stator. Different from the reverse flow in the rotor, which is 
caused by tip clearance vortex, the stator reverse flow is 
caused by suction surface flow separation. An animation of 
the time-varying solution (not shown here) shows that the 


stator separation zones fluctuate in size as the rotor passes by, 
but otherwise stay in the stator blades where they originated. 
The stator separation zones do not appear to rotate at this time. 
A similar observation was also noted by Davis and Yao 
(2006). 

At time E3, 1.1 revolutions after time E2, the flow formally 
enters the rotating stall. A zone of extended reverse flow cov- 
ering five rotor passages can now be seen in the front view of 
cut-plane 3 (330°) in figure 13(c) and (d). This zone sits in the 
middle of the less extensive reverse flow characterized by the 
green colored tip clearance flow on cut-plane 3. This less 
extensive reverse flow correlates with the 1 st mode. The 5- 
passage reverse flow zone initiates a spike-like disturbance on 
top of the 1 st mode and finally leads to the rotating stall. 

Further analysis of the simulation also indicates the correla- 
tion between the 1 st mode disturbance of axial momentum and 
the appearance of spike. The circumferential profile of the 
axial momentum near the rotor tip leading edge is low-pass- 
filtered to remove high-frequency noise so only the 1 st mode 
stays. Figure 14 shows the results at time El, E2, and E3. The 
magnitude of the 1 st mode disturbance grows from El to E3 as 
it rotates around the annulus. The locations of the trough of 
the waves (marked in figure 14) are where the flow has the 
lowest axial momentum. These locations match the regions of 
the larger reverse flows in all three cut-planes mentioned 
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Figure 14. — First mode disturbances of 
axial momentum at rotor tip. 
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(a) Stable condition (b) Near stall condition (time El) 
Figure 15. — Pressure contour near casing. 


previously. This is consistent with the observation of Day 
(1993) that spikes can grow out of the trough of the axial 
velocity modal waves if critical incidence is reached. 

As we pointed out before it is difficult to discern the actual 
transition from the rotating modal wave to the start of the 
rotating stall in figure 2. The simulation provides the details to 
complement the experiment. It is demonstrated here that a 
spike of 5 passages long does occur. It evolves smoothly out 
of the rotating mode and eventually leads to the breakdown 
(rotating stall). 


Interaction of Shock and Tip Clearance Flow 

Interaction between the tip clearance flow and shock in 
transonic compressor has been cited as one of the initiators of 
spike disturbance. As noted by Day (1993), short-length 
spikes can be viewed as embryonic stall cells with flow break- 
down in local regions; it is a 3D phenomenon. Computational 
simulations using a single blade passage to study the effect of 
the tip clearance flow on the onset of stall were reported by 
Adamczyk et al. (1993). One generally agreed conclusion is 
that the tip clearance vortex moves forward from stable to stall 
condition. Using vortex kinematics arguments, Hoying et al. 
(1999) developed a criterion for stall inception when the tra- 
jectory of the tip clearance vortex becomes perpendicular to 
the axial direction. Vo et al. (2005) also demonstrated that 
stall inception is most likely accompanied by the forward 
spillage of tip clearance vortex. One might ask: How does the 
shock interact with the tip clearance flow? Does the shock 
directly cause flow separation that may contribute to the stall 
inception? 

To answer these questions, we take a closer look at the flow 
within the rotor passage. A comparison of pressure contours 
near the casing between the stable and near-stall (time El) 
condition is made in figure 15. In the stable condition, the 
pressure fields in the two passages shown are almost identical. 
There exists a distinguishable attached oblique passage shock, 
which originates from the leading edge of the rotor blade and 
extends across the passage to 80% chord on the suction sur- 
face of the adjacent blade. At time El, the pressure fields 
between passages are no longer identical. Even though the 
flow structure is not identical among passages there is still 
similarity between them before the flow breakdown. We will 
discuss the common features that are generic for all passages 
before stall. Examining the pressure field at time El, we see 
the single attached oblique shock now changes to a more 
complicated dual-shock system: a strong detached shock fol- 
lowed by a weaker secondary shock. The strong detached 
shock, together with the stagnant flow around the leading 
edge, creates a high pressure region near the leading edge of 
the pressure surface. The pressure outside this area is rela- 
tively low until aft of the secondary shock where the pressure 
level once again rises but now fills uniformly between the 
blades. 

Streamlines forming the tip clearance vortex for both stable 
and near-stall (time El) conditions are shown in figure 16. The 
pressure contour near the casing is overlaid in the background to 
help the discussion. First take a look at the stable condition. 
Only the forward 15% portion of the tip clearance flow emerges 
to form the core of the tip clearance vortex. The trajectory of the 
tip clearance vortex is not perpendicular to the axial direction. 
The oblique shock intersecting the vortex only slightly alters the 
shape of the vortex. The vortex moves away from pressure 
surface of the adjacent blade without impinging it. 
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(a) Stable condition (b) Near stall condition (time El) 

Figure 16. — Stream lines forming the tip clearance vortex. 



(a) Stable condition 



(b) Near stall condition (time El) 

Figure 17. — Velocity vector on a meridional 
plane 30% pitch off pressure surface. 


At the near-stall condition (time El), figure 16(b), the core 
vortex, again formed by the forward 15% of the tip clearance 
flow, now moves forward. The center portion of the vortex 
trajectory becomes perpendicular to the axial direction. This 
trend follows well with the stall inception criterion of Hoying 
et al., (1998). However, the shape of the vortex trajectory is 
first pushed backward by the high pressure near the pressure 
surface leading edge and then forward by the high pressure 
behind the secondary shock as it proceeds downstream. This 
distortion of the vortex trajectory is caused by the dual-shock 
system. The criterion of Hoying, developed for low- speed 
machines, still holds for areas not affected by the dual-shock. 
It is unclear whether or not the distorted vortex can accelerate 
the flow breakdown, but it is clear that the distortion occurs as 
the flow is close to stall with the existence of the dual-shock. 

Flow is likely to separate due to the shock-boundary layer 
interaction, in which the low momentum flow in the boundary 
layer is subjected to the sudden pressure rise across the shock. 
Thus, it is of interest to examine whether shock- induced flow 
separations occur in the simulation. There are two possible 
locations for flow separation in the rotor passage: one is on the 
suction surface and the other is on the casing surface. By 
examining the streamlines in figure 1 6, we rule out the suction 
surface separation. To search for the casing separation, the 
velocity vector in the meridional view is examined. 
Figure 17(a) shows such a plane 30% pitch off the pressure 
surface. There is no observable reverse flow in this view. This 
indicates that the vortex moves straight downstream and that 
the flow remains attached to the casing. At time El, the re- 
verse flow associated with the distorted vortex can be seen in 
figure 17(b). The velocity vector shows that the flow reversal 
occurs beneath the casing and is caused by the organized spi- 
raling fluids of the distorted vortex. There is no indication of 
casing separation at this instant. The cross-section of the vor- 
tex, however, is larger than the stable condition. This also 
shows that the negative axial velocity of the rotor flow shown 
in figure 13 is caused by the spiraling fluids of the tip clear- 
ance vortex and has no association with casing or blade sepa- 
ration. With the above observations, we conclude there is no 
shock-induced separation in this during the stall inception. 

As the compressor proceeds more into stall between times 
E2 and E3, the flow is breaking down, i.e., the mass flow rate 
through the compressor is declining. The reverse flow in a few 
passages combines to form the large blockage as shown in 
figure 13(b) and (c), which in turn expels flow to adjacent 
blades. The blockage creates the spike-type disturbance over a 
few passages. The expelled flow entering the adjacent blade 
passages brings in extra axial momentum to the casing areas 
of these passages. The additional axial momentum behaves 
like tip injection, which reduces rotor incidence levels, relative 
to time El. The vortex shape and shock structure of two pas- 
sages at time E2 are compared in figure 18. Figure 18(a) is a 
passage where the spike occurs. The shape of the vortex tra- 
jectory and shock structure is similar to those at time El, 
figure 16(b), only that the vortex continues to grow and the 
high pressure behind the detached shock pushes the vortex 
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backward so the trajectory moves away from being normal to 
the axial direction. Spillage of tip clearance forward of the 
leading edge can be seen in figure 19(a), where it is colored by 
the pink stream-lines. In the passage where the flow condition 
improved, as shown in figure 18(b), the dual-shock system 
returns to a single- shock and the vortex trajectory is not dis- 
torted. There is no spillage of the tip clearance flow forward of 
the leading edge, figure 19(b). 




(a) Spillage in (b) No spillage in 

deteriorating passage improving passage 

Figure 19. — Tip clearance flow showing spillage 
over rotor leading edge at time E2. 


Simulation Summary 

(1) The simulation results of the high speed axial compres- 
sor Stage 35 provide evidence that the current numerical pro- 
cedure may be adequate for predicting the onset of flow 
instabilities and their subsequent growth into a fully developed 
rotating stall. The simulation results clearly show a full- 
annulus 1 st mode disturbance starting to travel at 100% rotor 
speed at the peak of the pressure rise transitions to a spike 
disturbance that traveling at 84% rotor speed. This spike dis- 
turbance then splits into multiple stall cells while their speeds 
reduce. Finally the multiple stall cells coalesce into a large 
single stall cell that rotates at 43% of rotor speed when fully 
developed. This simulation models the natural growth of the 
stall process without the use of an artificial inlet distortion. 

(2) As the simulation is throttled to stall, small amplitude 
long-length waves of 1 st mode and high-order modes between 
the 1 st and the 36 th can be discerned. The magnitudes of the 
high-order modes fluctuate, but do not grow much during the 
transition to stall. The 1 st mode wave, however, persistently 
grows before breakdown. A short-length disturbance covering 
five rotor passages bursts into a spike-like disturbance during 
the breakdown. This spike disturbance then grows into a rotat- 
ing stall. 

(3) The computed 1 st mode disturbance has a length scale 
closely resembling the less than one rotor circumference dis- 
turbance prior to stall in the experiment. The computed speed 
of the 1 st mode disturbance (100% rotor speed) is also similar 
to the measured speed (80% rotor speed) of the 1 st mode dis- 
turbance from the experiment, estimated from the results of 
figure 2. The difference is within the resolution of the picture 
quality of the published data. 

(4) The 1 st mode of the axial momentum correlates well 
with the appearance of the growing reverse flow. The trough 
of the axial momentum matches the region of the larger re- 
verse flow that leads to the birth of the initial spike. 

(5) The interaction of the shock and tip clearance vortex is 
examined. As the compressor approaches stall the trajectory of 
the tip clearance vortex moves forward to become perpendicu- 
lar to the axial direction. The shock system, in the meantime, 
also moves forward, changing from an attached oblique sin- 
gle-shock to a dual-shock. The pressure field of the dual-shock 
system distorts the trajectory of the tip clearance vortex, creat- 
ing a large reverse flow region in the distorted vortex. The 
reverse flow in the rotor is the outward tip clearance vortex as 
no shock- induced casing- or blade-separation is found. 

(6) During the flow breakdown, spillage of the tip clearance 
flow forward of the rotor leading edge is found in several 
passages, where the flow condition is deteriorating. A large 
blockage region is formed in the deteriorating passages; this in 
turn impels flow to enter into adjacent passages. The influx of 
fluid momentum in the adjacent passages improves the flow 
condition there. 
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Conclusions 

The TURBO code has been applied to the 3D full-annulus 
simulation of a transonic compressor stage NASA Stage 35 to 
investigate the flow evolution prior to stall. This simulation 
does not include the system volume, and therefore cannot 
model the surge mode oscillation that arises from the test 
facility ducting system. With the surge mode removed, the 
stalling flow first takes the form of a rotating long-length 
(modal) disturbance; it then changes into a spike. The spike 
seems to grow out of the long-length disturbance and eventu- 
ally leads to stall. The spike emerges when the tip clearance 
vortex becomes perpendicular to the axial direction and in 
some passages, with forward spillage. The passage shock 
changes from single-shock to dual-shock when approaching to 
stall. The dual-shock distorts the trajectory of the vortex but 
doesn’t cause flow separation that directly contribute to the 
stall. The simulation correlates well with the experiment data 
excluding the surge mode. The computational resources in- 
volved in conducting such a simulation are prohibitively large 
at the time and are not suitable for routine executions. While 
this is certainly an issue, the key point is that the simulations 
of this type generate qualitatively and quantitatively accurate 
flow field data which can support our understanding and infer- 
ences of the unsteady compressor flow field, including support 
and interpretation of experimental work. 
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